Bell Inequalities for Graph States 
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We investigate the non-local properties of graph states. To this aim, we derive a family of Bell 
inequalities which require three measurement settings for each party and are maximally violated 
by graph states. In turn, for each graph state there is an inequality maximally violated only by 
that state. We show that for certain types of graph states the violation of these inequalities in- 
creases exponentially with the number of qubits. We also discuss connections to other entanglement 
properties such as the positivity of the partial transpose or the geometric measure of entanglement. 
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Quantum theory predicts correlations which are 
stronger than the correlations of local hidden variable 
(LHV) models. By definition, LHV models have to obey 
the constraints of realism and locality: Any observable 
has a predetermined value, regardless of whether it is 
measured or not, and the choice of which observable to 
measure on one party of a multipartite system does not 
affect the results of the other parties. These constraints 
lead to the so-called Bell inequalities which put bounds 
on the correlations. These inequalities turn out to be 
violated by certain quantum mechanical states Q, ■ 

In this letter we address the question, whether graph 
states allow a LHV description or not. Graph states 
form a family of multi-qubit states which comprises many 
popular states such as the Greenberger-Horne-Zeilinger 
(GHZ) states and the cluster states |2j • Graph states are 
also crucial for applications: All codewords in the stan- 
dard quantum error correcting codes correspond to graph 
states and one-way quantum computation uses graph 
states as resources Recently, graph states have been 
produced in optical lattices p| and the basic elements 
of one-way quantum computing have been demonstrated 
experimentally 0- Also, general methods for the gener- 
ation of graph states have been explored @ . 

It is a natural and important question whether these 
tasks and experiments, including the effects of noise can 
be described by LHV models. To answer this question, 
we derive a class of Bell inequalities. Each graph state vi- 
olates one of these inequalities in the GHZ sense, i.e., by 
saturating all correlation terms and for certain types of 
graph states the violation of local realism increases expo- 
nentially with the number of qubits. In this way we show 
that tasks like measurement based quantum computing 
and quantum error correction are far from the realm of 
LHV theories. Note that the non-locality of special ex- 
amples of graph states has been shown recently . 

Graph states are defined as follows. Let G be a 



FIG. 1: Types of graphs for the case of five vertices: (a) The 
linear cluster graph LC5. (b) The ring cluster graph RC5. (c) 
The star (or GHZ) graph ST 5 . This describes a GHZ state, (d) 
The fully connected graph FC5. This graph can be obtained 
from ST5 by local complementation on the second qubit (see 
Lemma 4). It also describes a GHZ state 0. 



graph, i.e., a set of n vertices and some edges connecting 
them. Some interesting graphs are shown in Fig. 1. For 
each vertex i the neighborhood N(i) denotes the vertices 
which are connected with i. We can associate to each 
vertex i a stabilizing operator gi by 
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From now on, XM,V«,ZW denote the Pauli matrices 
cr x ,a y ,a z , acting on the i-th qubit. For instance, for 
the fully connected three vertex graph, the stabilizing 
operators are g x = X^ Z^ Z^\ g 2 = Z^X^Z® and 
g 3 = ZWZWXW. The graph state \G) associated with 
the graph G is the unique n-qubit state fulfilling 



9i \G) = \G), fari = l, 



(2) 



Physically, the graph describes the perfect correlations 
in the state \G), since (g l ) = (X® ® jeN{l) = 1. At 
the same time, it denotes a possible interaction history 
leading to \G), i.e., \G) can be produced by an Ising type 
interaction acting between the connected qubits. 

Given the stabilizing operators gi, we can look at the 
group of their products, the so-called stabilizer 



5(G) = { Sj ,j = 1,...,2"}; 
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where Ij(G) denotes a subset of the vertices of G. If a 
certain generator g^ appears in the product in Eq. J3J, 
(i.e., k € Ij{G)), we say that Sj contains g^. The group 
S(G) is commutative and has 2™ elements. Of course, for 
all these elements Sj\G) = |G) holds. A last property we 
have to mention is the fact that 



i=l 



2 n \G)(G\ 



(4) 



holds, as can be checked by direct calculation 

Now we present our idea for the derivation of Bell in- 
equalities. Given a graph G all the stabilizing operators 
are of the form 
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where the single qubit observables are either the 
identity or one of the Pauli matrices: O 
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{l,±X( fc \±y( fc ),±Z( fc )}. We look at the operator 

2™ n 



B(G) = J2MG) = Y, < S>° i l 
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To give a simple example, this operator for the fully con- 
nected graph state for three qubits reads 

B(FC 3 ) = l^l^l^ +X^Z^Z^ +Z^X^Z^ 
+ Z {1) Z {2) X {3) + Y {1) Y {2) 1 {3) + y( 1 )i( 2 )y< 3 ) 
+ i<l)y(2)y(3) -X (1) Z (2) X (3) . (7) 

We take this 2? as the Bell operator and compute a bound 



C(G) = C(B) = max|{B}|, 

LHV 



(8) 



where the maximum of the absolute value of the mean 
value (B) is taken over all LHV models. Here, it suf- 
fices to look at deterministic LHV models which have to 
assign definite values {+1,-1} to the observables 

whenever 0$ ^ 1 @ ■ This is due to the fact that nonde- 
terministic LHV models can be viewed as deterministic 
LHV models where the hidden variables are not known. 
In principle, in the definition of B the Pauli matrices can 
be replaced by arbitrary dichotomic observables. Since 
we are interested in graph states, we will, however, al- 
ways use X, Y, Z. 

If we can find for a given graph G a bound C(G) < 2™, 
the nonlocality of the graph state \G) is detected. This 
is due to Eq. 10}, which implies that for the graph state 
(B) = 2" holds. Also, the graph state violates the Bell 
inequality maximally. In the example of Eq. {J} we will 
see later that C(FC 3 ) = 6. This gives rise to the Bell 
inequality \(B(FC3))\ < 6 which is violated by the state 
I-FC3). In the following, it will also be useful to compare 



the strength of the Bell inequalities by the normalized 
parameter V{G) := C{G)/2 n or 1/V. Note that we have 
a valid Bell inequality whenever T> < 1 . 

So the main task is to find the value of C(G) or 
T>(G). An exact calculation is, in general, very demand- 
ing. However, as we will show, it is quite easy to obtain 
bounds on T>(G) depending on the structure of the graph 
G, especially when we can identify some subgraphs where 
the bounds are already known. Together with the exact 
calculation of T>(G) for graphs with a small number of 
qubits this allows us to derive some general results for 
arbitrary graphs. Let us first note two useful facts about 
the dependence of V on the LHV models. 

Lemma 1. We can restrict our attention to LHV mod- 
els which assign +1 to all Z measurements. 
Proof. In an element Sj of the stabilizer we have £ 
{yW, Z«} iff the number of Y and X in N{i) is odd. So 
if a LHV model assigns —1 to Z^ l \ we can, by changing 
the signs for Z^,Y^ and for all X^ and all Y^ with 
k € N(i), obtain a LHV model with the same mean value 
of B and the desired property. □ 

Lemma 2. Let B be a Bell operator for an arbitrary 
graph and let B' be the Bell operator which is obtained 
from B by making a permutation 

P : {lW,JCW,yW,2W} ^ {i(i) iX ^,Y^,Z^} (9) 

of the observables on one qubit. Then T>(B) = V(B'). 
Proof. It suffices to show the Lemma for a transposi- 
tion OW <_> OW, Some transpositions are a renaming 
of variables, and the interesting transpositions are of 
the type AW <-> lW with A® = jW, yW, say 

= A« for definitcness. If a given LHV model 
LHVi assigns +1 to X^ the transposition X^ <-> 1^ 
does not change (B). If the LHV model assigns —1 to 

we can construct a new LHV model LHV2 from 
LHVi by flipping the signs from Y^,Z^. This fulfills 
{B')lhVi — —{B)lhv 2 - This proves the claim, since C is 
defined via the absolute value. □ 
Now we derive an estimate for T>(G), when G is built 
out of two other graphs G\ and G2 in a certain way. 

Lemma 3. Let G\,Gi be two graphs and let G be 
the graph which comprises G\ and G2 and one single 
connection between one vertex of G\ and one of G2, i.e., 
G = Gi — G 2 . Then 



V(G) <V(G 1 )V{G 2 ). 



(10) 



Proof. The proof is given in the Appendix. □ 
It is much more demanding to derive bounds on T>(G) 
when G is made out of subgraphs in a more compli- 
cated way than the way above. However, it is easy to 
see that T>(G) < 1 whenever G contains a subgraph G\ 
with T>(Gi) < 1. This is due to the fact that the stabi- 
lizer of Gi is a subset of the stabilizer of G up to some 
extra Z terms which can be neglected due to Lemma 1. 

Finally, we want to show the invariance of T> under 
the so-called local complementation of a graph. This 
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TABLE I: Value for T>{G) for different interesting graphs (see 
Fig. 1) and different numbers of qubits n. 



n 


3 


4 


5 


6 


7 


8 


9 


10 


LCn 


3/4 


3/4 


5/8 


9/16 


8/16 


7/16 


25/64 


22/64 


RC„ 


3/4 


3/4 


5/8 


7/16 


7/16 


6/16 


21/64 


19/64 


ST n 


3/4 


3/4 


5/8 


10/16 


9/16 


9/16 


34/64 


34/64 


FC n 


3/4 


3/4 


5/8 


10/16 


9/16 


9/16 


34/64 


34/64 



transformation acts as follows: One picks up a vertex 
io and inverts the neighborhood N(io), i.e., all connec- 
tions between two vertices belonging to N(iq) are cut and 
vertices in N(io) which were unconnected become con- 
nected. Connections between N(io) and the rest of the 
graph are not affected. To give an example the graph 
ST n can be transformed by a local complementation on 
the central qubit into the graph FC n (see Fig. 1). The 
local complementation of a graph acts on the graph state 
as a local unitary transformation of the (local) Clifford 
group 0, This means that it transforms on each 

qubit Pauli matrices into Pauli matrices. So we have: 

Lemma 4. Let G\ be a graph and G-i be a graph 
which arises from G\ by local complementation. Then 
V{G 1 )=V{G 2 ). 

Proof. Since the local complementation maps Pauli ma- 
trices to Pauli matrices on each qubit, T> is not changed 
due to Lemma 2. □ 

After showing in the previous Lemmata how to esti- 
mate T>(G) when the values of subgraphs are given, it 
is now time to calculate T>(G) for graphs with a small 
number of qubits n. This can easily be done by a com- 
puter, by checking (£>) for all the 8™ LHV models. In 
these calculations, Lemma 1 reduces the computational 
complexity significantly. The numerical results for inter- 
esting graphs up to ten qubits are shown in Table I. 

Let us shortly discuss these results. First, note that 
all of the states in the table violate a Bell inequality of 
the type |(B)| < C{G) since for all states and n in Table I 
T>(G) < 1 holds. Then, it is remarkable, that long linear 
chains and large rings have a small T> (i.e. the violation 
of the Bell inequality is large), while the violation for the 
GHZ state is not so large. This in contrast to the usual 
Mermin inequality for several qubits. Due to Lemma 4, 
the values V{ST n ) and V(FC n ) always coincide. How- 
ever, it is interesting that V(FC 6 ) £ V(FC 3 )V(FC 3 ), 
thus a generalization of Lemma 3 to arbitrary connec- 
tions between the graphs is not true. In general, we have: 

Theorem 1. Any graph state violates local realism. 
Proof. If the graph consists only of two connected ver- 
tices, the graph state is equivalent to a two qubit sin- 
glet state, which violates the original Bell inequality 
Connected graphs with more vertices always contain a 
subgraph with three vertices. Due to Table 1 and the ar- 
gumentation after Lemma 3 this implies that the graph 
state violates local realism. □ 



Theorem 2. Let G±, G2, G3, ■■■ be a family of tree 
graphs (i.e., graphs which do not contain any closed ring) 
with an increasing number of vertices such that each Gi 
contains as a subgraph a linear chain of a size which in- 
creases linearly with i. Then violation of a Bell inequality 
for \Gi) increases exponentially with i. 
Proof. From Lemma 3 and Table I it follows that for 
linear cluster graphs LCi the value l/T>(LCi) increases 
exponentially. For tree graphs Lemma 3 can be applied 
again to show that T>(Gi) is smaller than the T> from the 
longest linear cluster in Gi. □ 

Similar statements as the one expressed in Theorem 2 
can be derived for other families of graphs. In fact, all 
graphs which are built by connecting an increasing num- 
ber of subgraphs via single edges exhibit an exponentially 
increasing violation of local realism, due to Lemma 3. 

Let us compare the Bell inequalities with other en- 
tanglement properties. Here, our restriction to the ob- 
servables X, Y, Z becomes crucial. A Bell inequality 
I < C(G) is equivalent to a witness of the type 

W = V(G)-\G)(G\, (11) 

i.e., a quantum mechanical state violates the Bell in- 
equality iff (W) < This has two conse- 
quences. First, since all fully separable states = 
\a)\b}...\n) admit a LHV description, here T>(G) > 
max| a \|M . .| n ) I ( a i b, n |G) | 2 has to hold. The quantity 
1 — maxi \|(,\ i n \ \{a,b, ...,n\G)\ 2 has been shown to be 
an entanglement monotone for multipartite systems, the 
so-called geometric measure of entanglement [12| . So our 
bounds on C also deliver lower bounds for this measure of 
entanglement for graph states. In turn, the fact that the 
geometric measure equals 1/2 for all GHZ states implies 
that always V(ST n ) > 1/2. Second, we can state: 

Theorem 3. If V(G) > 1/2 then the Bell inequality 
|(B) I < C(G) detects only states which have a negative 
partial transpose with respect to each partition. 
Proof. Let us fix a bipartite splitting for the multipartite 
system. The graph state has a Schmidt decomposition 
\G) = X)i=i a iN) w hh respect to this splitting. It is 
known that then 1/fc < a 2 , < 1/2 if ao is the biggest 
Schmidt coefficient 13]. If we define \ip) = \H)/Vk 
it is also known that the witness VV = 1/k — |V')('0| de- 
tects only states which have a non-positive partial trans- 
pose with respect to this partition (lif. However the 
witness W detects even less states, since W — fcao 2 W = 
(V-al)l-\G){G\ + kal\i>)^\>kal\jf)^\-\G){G\>Q 
which implies that (W) < only if (VV) < 0. □ 

In conclusion, we have derived a family of Bell inequal- 
ities for multipartite systems based on the graph state 
formalism. These inequalities are maximally violated by 
graph states and allow to detect the non-locality of all 
graph states. Also, the inequalities can be related to 
other topics as the geometric measure of entanglement 
and the criterion of the partial transposition. The fact 
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that graph states do not admit a LHV model strongly 
suggests that tasks like measurement based quantum 
computation and quantum error correction cannot be de- 
scribed within classical physics. 

We thank Marc Hein for discussions and P.H. thanks 
the Universitat Innsbruck for hospitality. This work has 
been supported by the DFG (Schwerpunkt Quantenin- 
formationsverarbeitung) and the EU (MEIF-CT-2003- 
500183, IST-2001-38877, -39227, OLAQI). 

Appendix — Here we give the proof of Lemma 3. Let 
us consider a given LHV model and assume that the con- 
nection is between the vertices «o € Gi and jo G G 2 . We 
can write any stabilizing operator of G as sy(G) = aibj 
where contains only g\~ with k S G\ and bj con- 
tains only gk belonging to G2. We can arrange the set 
{cii} — m U ... U ci6 into six disjoint subsets in the follow- 
ing way: Oi, 02 contain the a,i where gi is absent, 03, 04 
have an X^ ' at the vertex io and 05,06 have a Y^°\ 
Note that this implies that 03, 04, as, and ci6 comprise the 
di which contain gi . For dj € 01,03,05 the LHV model 
gives (di) = +1 while for the a, G 02,04 or 06 we have 
(di) — — 1. Let us denote the number of elements in the 
sets Oi, o 2 , o 3 , 04, o 5 and o 6 by p a , q a , r a ,s a ,t a and u a . We 
can define j a =: p a - q a , S a := r a - s a and e a = t a - u a . 
Then we have \{B(G{))\ = | 7o + S a + e a \ < C(Gi). We 
can do the same ordering for the bj, introducing bi, b6 
and p b , ...,Ub and 76) S b ,s b , with \j b + S b + e b \ < C(G 2 ). 

Let us investigate the relationships between 7, 5, and 
e in some more detail. By flipping the sign which is 
assigned to V^ ' by the given LHV model we can con- 
struct a new LHV model with | (B(G\))\ = \j a — fa+£ a \ < 
C(Gi). We can also flip the signs of Y^°\X^ or F<» 
leading to new bounds of the type I7 ± 5 ± e\ < C. All 
in all, this leads to the conclusion, that for all sixteen 
combinations of signs 

|( 7a ±S a ± e a )( 7b ±S b ± e b )\ < C{G 1 )C(G 2 ) (12) 



signs. They can be written in the following way: 



holds. Finally, note that the operators Sij can be grouped 
into 36 groups according to s%j = <Xibj. The mean values 
(sij) = (aibj) would factorize if there were no connection 
between the graphs. In this case, the Lemma is trivial. 

What changes for the (s^) due to the extra con- 
nection? The (sij) can be written in a 6 x 6 block 
matrix according to the grouping into the Sy, where 
each block bears the sign of the corresponding (sy). 
In the blocks Sij with i < 2 or j < 2 the extra 
connection only introduces transformations of the type 
Z <-> 1 at io or jo, which can be neglected due to 
Lemma 1. More interesting is the 4x4 block matrix 
of the blocks s y - with 3 < i,j < 6. One can calculate 
that here the extra connection induces the transforma- 
tion {X^X^ \X^Y^ \Y^X^°\Y^Y^} 1 — > 
{Y( l °)Y { i°\-Y {i ^ X { 3°\-X^Y { i°\ X^XUo)} on the 
vertices io and jo- So, depending on the LHV model, this 
results in this block matrix in two possible changes of 
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— h 
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+ - 
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. — h 


h . 



" -1 

- + 
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+ - 


h 

. H 


+ - 
h . 
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For the first possible transformation we have to show 
that I J2ij( s ij)\ = \la{lb + S b + e b ) + <5 Q (7h + 5 b - e b ) + 
£ a {jb — $b +£&)| < C(Gi)C(G 2 ) This can be derived from 
Eq. (|T2*|l . distinguishing 64 cases depending on the signs 
of j a ,Jb,S ai S b ,e a and e b . For instance, if j a ,S a ,S b > 
and 76, e a , e b < we use (70 + S a + e a )("/b + 5b + e&) < 

7a(7b + Sb + £b) + 5 a {"fb + $b - Eft) + £a {"fb ~ fa + Efe) = 

J2ij( s ij) < -( 7a + (5 a -e a )(7fc-(5 fc + e fc ), yielding an upper 
and a lower bound for J2ij( s ij)- The proof of the other 
63 cases and the second transformation is similar. □ 
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